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SUMMARY 



Pluripotent stem cells (PSCs) occupy a spectrum of reversible molecular states ranging from a naive ground-state in 2i, to metastable em- 
bryonic stem cells (ESCs) in serum, to lineage-primed epiblast stem cells (EpiSCs). To investigate the role of DNA methylation (5mC) 
across distinct pluripotent states, we mapped genome-wide 5mC and 5-hydroxymethycytosine (5hmC) in multiple PSCs. Ground-state 
ESCs exhibit an altered distribution of 5mC and 5hmC at regulatory elements and dramatically lower absolute levels relative to ESCs in 
serum. By contrast, EpiSCs exhibit increased promoter 5mC coupled with reduced 5hmC, which contributes to their developmental 
restriction. Switch to 2i triggers rapid onset of both the ground-state gene expression program and global DNA demethylation. Mecha- 
nistically, repression of de novo methylases by PRDM14 drives DNA demethylation at slow kinetics, whereas TETl/TET2-mediated 5hmC 
conversion enhances both the rate and extent of hypomethylation. These processes thus act synergistically during transition to ground- 
state pluripotency to promote a robust hypomethylated state. 



INTRODUCTION 

Pluripotent stem cells (PSCs) exhibit the core properties of 
self-renewal and the capacity to differentiate into all 
embryonic lineages. These properties are underpinned by 
expression of a core circuitry of pluripotency transcription 
factors and reinforced by epigenetic mechanisms (Surani 
et al. ( 2007). Mouse embryonic stem cells (ESCs) are classi- 
cally cultured in serum/leukemia inhibitory factor (LIF) 
and considered to harbor "naive" pluripotency, which 
reflects their unrestricted potential to contribute to all 
somatic lineages and the germline in chimeras (Nichols 
and Smith, 2009). However, the undefined signaling in 
serum conditions means that ESCs exhibit heterogeneous 
gene expression profiles and exist in at least two distinct 
but interchangeable populations corresponding to a naive 
and a more developmentally restricted "primed" state and 
are thus termed metastable (Hayashi et al., 2008; Toyooka 
et al., 2008). 

Recent studies have revealed that inhibition of mitogen- 
activated protein kinase kinase (MEK) and glycogen syn- 
thase kinase-3 (GSK3) by culture in 2i harnesses ESCs in 
an apparently homogeneous naive "ground-state," which 
more closely recapitulates the transcriptional and epige- 
netic state of the preimplantation epiblast cells from which 
ESCs are derived (Ying et al., 2008; Marks et al., 2012; Leitch 
et al., 2013). In contrast, epiblast stem cells (EpiSCs), which 
are established from the postimplantation epiblast, occupy 
a pluripotent state "primed" toward lineage differentiation, 



which is reinforced by culture with FGF and activin (Brons 
et al., 2007; Tesar et al., 2007). An additional PSC type 
termed embryonic germ cells (EGCs) can be derived from 
primordial germ cells (PGCs). Although PGCs are specified 
from primed postimplantation epiblast, EGCs appear to be 
transcriptionally and functionally comparable to ESCs and 
not to EpiSCs, implying that EGCs undergo reprogram- 
ming to a more naive state (Matsui et al., 1992; Leitch 
et al., 2010, 2013). Collectively, these mouse PSCs repre- 
sent a wide spectrum of pluripotent states that can be 
used for investigation of developmental processes 
in vitro. Moreover, understanding the molecular determi- 
nants of naive pluripotency in mouse ESCs may yield in- 
sights toward derivation and benchmarking of bona fide 
naive human ESCs (hESCs) (De Los Angeles et al., 2012). 

Epigenetic mechanisms including histone modifications 
and DNA methylation play an important role in stabilizing 
cell identity and orchestrating many developmental pro- 
cesses. DNA methylation (5mC) is associated with tran- 
scriptional silencing and is essential for mouse embryonic 
development but dispensable for the derivation of ESCs 
(Okano et al., 1999). Nonetheless, DNA methylation-defi- 
cient ESCs are unable to differentiate, indicating that 
5mC is required for exit from pluripotency (Jackson et al., 
2004). Indeed, de novo promoter methylation occurs rela- 
tively early during differentiation — at transition to progen- 
itor cells rather than upon terminal differentiation — which 
may reflect an early and stable "lock" on lineage commit- 
ment (Mohn et al., 2008). Moreover, 5mC at key gatekeeper 
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genes such as Elf5 acts to restrict ESCs to embryonic fates 
and prevent their contribution to extraembryonic tissues 
(Ng et al., 2008). Notably, metastable ESCs in serum exhibit 
relatively high global DNA methylation levels that are 
comparable with differentiated tissues, whereas naive 
ground-state ESCs exhibit a globally hypomethylated state 
(Leitch et al., 2013; Yamaji et al., 2013; Ficz et al., 2013; 
Habibi et al., 2013). Thus, the genomic landscape of DNA 
methylation may play an important role in denning the 
balance between restriction and plasticity in PSCs. 

Recent studies have shown that 5mC can be oxidized to 
5-hydroxymethycytosine (5hmC) by TET enzymes and 
that this process may play a role in global DNA demethyla- 
tion during preimplantation development and in the germ- 
line (Tahiliani et al., 2009; Iqbal et al., 2011; Hackett et al., 
2013). Conversion of 5mC to 5hmC can potentially 
engender both replication-dependent "passive" or replica- 
tion-independent "active" demethylation (Hackett et al., 

2012) . These pathways could allow turnover or resetting of 
DNA methylation to promote and maintain pluripotent 
states. Indeed, both 5hmC and TET enzymes are relatively 
abundant in PSCs (Ficz et al., 2011; Koh et al., 2011). How- 
ever, Tetl/Tet2 null ESCs can be efficiently derived and are 
pluripotent (Dawlaty et al., 2013), and thus, any potential 
role of 5hmC in promoting dynamic DNA methylation 
turnover in the context of each pluripotent state remains 
unclear. Indeed, alternative DNA demethylation mecha- 
nisms involving suppression of de novo or maintenance 
methylation activity have been suggested (Yamaji et al., 
2013; Seisenberger et al., 2012), implying that multiple par- 
allel mechanisms may promote DNA demethylation in PSCs 
(Hackett et al., 2012). Importantly, whereas recent studies 
have mapped the hypomethylated epigenome of ground- 
state ESCs at base resolution (Habibi et al., 2013; Ficz et al., 

2013) , they did not distinguish between 5mC and 5hmC. 
Moreover, the global profiles of 5mC and 5hmC in the 
primed state (EpiSC) and in EGCs remain to be determined, 
as do the precise mechanisms that drive interconversion be- 
tween distinct DNA methylation landscapes. 

Here, we set out to investigate genome-wide 5mC 
and 5hmC in diverse PSCs and found significant changes 
associated with distinct pluripotent states. Moreover, we 
show that dual TET-mediated 5hmC conversion and 
PRDM14-dependent Dnmt3 repression synergistically 
drive DNA demethylation during transition to ground- 
state pluripotency. 

RESULTS 

Global 5mC and 5hmC in Pluripotent Cells 

To investigate the roles and dynamics of 5mC and 5hmC 
in PSCs derived from distinct developmental origins, or un- 



der different culture conditions, we used live genetically 
matched cell types: male XYand female XX ESCs on feeders 
in classical serum/LIF media, ground-state ESCs and EGCs 
under serum-free 2i/LIF conditions (both XY), and primed 
EpiSCs derived from postimplantation epiblast. We first 
determined the global levels of 5mC and 5hmC quantita- 
tively by ELISA and found that total 5mC levels were highly 
comparable among XY ESCs in serum, EpiSCs, and control 
mouse embryonic fibroblasts (MEFs). However, 5mC levels 
were reduced slightly in XX ESCs and significantly reduced 
by 2.2-fold in both ESCs and EGCs maintained in 2i (Fig- 
ure 1A). The global levels of 5hmC followed a similar 
pattern, with relatively high levels in ESCs in serum and 
EpiSCs but reduced levels in cells in 2i. The global reduc- 
tion of 5hmC in 2i was comparable in extent to 5mC 
(1.9-fold); however, 5hmC is significantly higher in all 
PSCs relative to MEFs (Figure 1A). Thus, global 5mC and 
5hmC levels are comparable in ESCs in serum and primed 
EpiSCs but significantly and proportionately reduced in 
ground-state PSCs cultured in 2i, which is consistent with 
recent studies (Leitch et al., 2013; Ficz et al., 2013; Habibi 
et al., 2013). 

To confirm changes in DNA modification, we used 
bisulfite sequencing to assay repetitive and single-copy 
genomic regions. This revealed that LINE-1 elements, 
CpG-dense promoters at Dazl and Rhox9, CpG-poor 
promoters at Esrrb and Elf 5, an intronic region at Cull, 
and exonic sequences of Rexl are all significantly hypo- 
methylated in ESCs/EGCs in 2i relative to ESCs in serum 
or EpiSCs, suggesting that hypomethylation occurs at 
most genomic landmarks in 2i (Figure IB; Figure SI 
available online). A notable exception is IAP elements, 
which showed no significant difference regardless of the 
culture condition (Figure IB), a result consistent with 
IAPs being resistant to DNA hypomethylation in vivo 
(Hajkova et al., 2002; Hackett et al., 2013). However, IAPs 
did exhibit significantly reduced DNA methylation in XX 
ESCs, which is apparently a specific effect because tested 
promoter regions and LINE-1 elements were not differ- 
entially methylated between our low-passage XX and XY 
ESCs, and global levels only showed a modest reduction 
in XX ESCs (Figures 1A and IB). Interestingly, we observed 
that XX, but not XY, ESCs exhibited additional global 
hypomethylation coincident with transition to feeder- 
free culture, while maintaining XX status, implying that 
female ESCs may be more susceptible to DNA demethyla- 
tion as a consequence of in vitro culture parameters than 
male ESCs (Figure SI). In particular, feeder-free culture 
may promote global hypomethylation in XX ESCs (Habibi 
et al., 2013; Zvetkova et al., 2005). Our data suggest that 
PSCs acquire distinct DNA methylation profiles that reflect 
their developmental origin and response to culture 
conditions. 
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Figure 1. 5mC and 5hmC Distributions in Distinct Pluripotent States 

(A) Global levels of 5mC (blue) and 5hmC (red) in EGCs in 2i (EGC-2i), ESCs in 2i (ESC-2i), male ESCs in serum (ESC-M), female ESCs in serum 
(ESC-F), EpiSCs, and MEFs, as a percentage of total nucleotides. 

(B) Bisulfite sequencing of repetitive and single-copy loci in distinct PSCs. Open circles represent unmethylated CpGs, and closed circles 
represent methylated CpGs. 

(C) Hierarchical clustering of PSCs based on genome-wide 5mC and 5hmC distribution. 

(D) Correlation coefficients of 5mC and 5hmC between PSCs. 

(E) ProfiLe shows the distribution of 5mC and 5hmC over an averaged metagene in each pluripotent state. 

(F and G) Profiles of 5mC and 5hmC in PSCs at (F) Esrrb and (G) Dazl. Shaded area indicates promoter region. 
**p < 0.01 of biological triplicate experiments. 



The Distribution of 5mC and 5hmC in Pluripotent 
Cells 

To further investigate the DNA modification landscape 
in pluripotent cells, we performed genome-wide methyl- 
ated DNA immunoprecipitation sequencing (meDIP-seq) 
and hydroxymethylated DNA immunoprecipitation 
sequencing (hmeDIP-seq). To facilitate comparison with 
(h)meDIP profiles from low-number PGCs, which share 
certain attributes of pluripotency and undergo dynamic 
DNA demethylation, we used a low-input protocol with a 
high-fragment resolution of ~150 bp (Hackett et al., 
2013). Hierarchical clustering of DNA modification profiles 
revealed that pluripotent cells cluster into distinct groups 
that correspond to ground-state, metastable, or primed 
cells (Figure 1C). Indeed, we observed similar clustering 



based only on 5mC and 5hmC profiles at binding sites of 
key chromatin factors (Figure SI). Consistently, the highest 
correlation coefficients were observed between cells within 
the same culture condition, regardless of developmental 
source, such that ESCs and EGCs derived from the ICM 
and PGCs, respectively, were highly comparable in 2i 
conditions. In contrast, the greatest difference was seen 
between ESCs/EGCs in 2i and EpiSCs, which also corre- 
sponds to the biggest difference in developmental poten- 
tial (Figure ID). Consistent with global hypomethylation 
under 2i conditions, ESCs in serum were more correlated 
with EpiSCs than with ESCs in 2i. 

Analysis of 5mC enrichment across a metagene revealed 
a comparable distribution in each pluripotent state, 
whereby 5mC is enriched within genie regions and 
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strongly depleted from transcriptional start sites (TSSs) (Fig- 
ure IE). We observed that 5hmC is also depleted from the 
TSS but found significant peaks of 5hmC enrichment in 
the flanking proximal promoter regions in ESCs in serum 
and EpiSCs, which was absent in ground-state cells (Fig- 
ure IE). We did note that 5hmC was enriched at the TSS 
specifically at low CpG density promoters but only when 
also bound by TET1 (Figure SI). Notably, some studies pre- 
viously reported that 5hmC is enriched directly at the TSS 
in general but used a 300-500 bp hmeDIP resolution rather 
than the 150 bp used here, which may have enabled some 
"bridging" of the TSS by fragments that also overlapped 
flanking promoter/exonic regions enriched in 5hmC (Fig- 
ure S2). Examination of specific loci such as the Esrrb 
promoter showed that it is highly enriched with 5hmC in 
ESCs in serum (XY and XX), as previously reported by 
Ficz et al. (2011), and that this is lost in 2i, coincident 
with reduction of 5mC (Figure IF). In contrast, 5hmC 
enrichment was maintained at the Dazl promoter in 
ground-state cells, which may be due to the retention of 
some 5mC signal under 2i conditions. However, there 
was a complete loss of 5hmC in primed EpiSCs together 
with a concomitant increase in 5mC enrichment 
(Figure 1G). 

Differential Promoter Methylation 

To investigate the genome-wide changes between PSCs, we 
examined promoters for differential methylation using XY 
ESCs in serum as a reference. A previous study has shown 
that feeder-free XX ESCs are globally hypomethylated 
(Zvetkova et al., 2005). However, consistent with our obser- 
vation of only a modest global depletion (~10%) of 5mC in 
XX ESCs on feeders, we found only 69 high-confidence (p < 
0.01) differential 5mC promoters between XX and XY 
ESCs. Surprisingly, there were more promoters methylated 
in XX ESCs (n = 55) relative to XY ESCs (n = 14), albeit there 
were very few differentially methylated promoters (DMPs) 
in absolute terms, which suggests that male and female 
ESCs share a highly similar promoter methylome, at least 
when cultured on feeders (Figure 2A). Comparison of 
ESCs in serum with ESCs in 2i revealed 228 high-confi- 
dence DMPs, with most being hypomethylated in 2i condi- 
tions (n = 179). Gene ontology (GO) analysis indicated that 
promoters hypomethylated in ground-state 2i relative to 
serum are not associated with specific processes/pathways 
after Bonferroni correction, suggesting that promoter hy- 
pomethylation is a general rather than a targeted effect in 
2i. Similarly, the few promoters that apparently exhibited 
higher 5mC in 2i were not enriched for GO categories but 
did apparently include Lin28b, which can promote both 
self-renewal and differentiation pathways depending on 
context (Shyh-Chang and Daley, 2013). Relatively elevated 
DNA methylation at Lin28b is consistent with lower Lin28b 



expression in 2i (Figure S3) (Marks et al., 2012). Interest- 
ingly, rare single-copy loci that escape DNA demethylation 
in PGCs, such as Sfil and Srrm2 (Hackett et al., 2013), are 
also resistant to demethylation in 2i conditions, although 
Vmn2r29 is hypomethylated (Figure S3). 

Examination of promoter 5mC in EpiSCs revealed a strik- 
ing number of promoters (n = 668) that become signifi- 
cantly hypermethylated relative to ESCs in serum, with 
only 19 exhibiting lower 5mC (Figure 2A). Promoters that 
become methylated in EpiSCs are highly enriched for being 
bivalently or H3K27me3 marked in ESCs and include many 
pluripotency and germline-specific genes such as Stella, 
Rexl, and Tell (Figures 2B and 2C). In addition to acquisi- 
tion of promoter 5mC, these genes often acquire DNA 
methylation throughout their entire transcription unit in 
EpiSCs, which may be related to stable heterochromatini- 
zation of these loci (Figure 2C). In general, when collated, 
the DMPs between all PSCs are enriched for an intermedi- 
ate CpG density and depleted of CpG-poor promoters 
(Figure 2B). 

Taken together, this analysis suggests that there is a 
general progressive gain of promoter 5mC coincident 
with transition from ground state > metastable > primed 
states. Notably, the most significant acquisition of pro- 
moter 5mC occurs in primed EpiSCs, which is in contrast 
to global 5mC, where significant gain occurs between 
ground-state ESCs in 2i and ESCs in serum. Thus, in gen- 
eral, promoters in ESCs cultured in serum are resistant to 
acquiring 5mC despite a high global DNA methylation 
context. In contrast, a significant number of promoters 
acquire DNA methylation in EpiSCs, coincident with a 
constraint on cellular plasticity. 

Analysis of differentially hydroxymethylated promoters 
(DhMPs) showed that both ground-state cells and primed 
EpiSCs were associated with a significantly reduced num- 
ber of 5hmC -enriched promoters relative to ESCs in serum 
(Figure 2A). The observation that the highest relative num- 
ber of 5hmC-enriched promoters occurs in serum ESCs is in 
contrast to 5mC, where in general, enrichment of methyl- 
ated promoters follows a linear pattern that correlates with 
increasingly restricted cells (ground state > metastable > 
primed). The enhanced promoter 5hmC in ESCs in serum 
may contribute to the relative resistance of promoters to 
acquiring 5mC, despite the high global 5mC levels. Consis- 
tent with this, tracking 5mC and 5hmC at all promoters 
across pluripotent states showed that promoter 5mC 
enrichment peaks in EpiSCs, whereas 5hmC enrichment 
peaks in ESCs in serum, with a similar result at biologically 
defined CpG islands (Figure 2D) (Illingworth et al., 2010). 
This effect was particularly prominent at pluripotency 
network genes, which presumably need to remain methyl- 
ation-free in ESCs to ensure robust expression and there- 
fore may be enriched in 5hmC (Figure 2E). In contrast, 
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Figure 2. Differentially DNA-Modified Promoters between PSCs 

(A) Number of DMPs and DhMPs in PSCs relative to maLe ESCs in serum. Black bars indicate the number of promoters with higher DNA 
modification, and gray bars indicate the lower levels in the indicated cell type. 

(B) Genetic and epigenetic features of differentially methylated promoters (DMP). Upper panel shows CpG density and chromatin state of 
DMPs that exhibit gain of 5mC between each successive pluripotent state (2i to serum, serum to EpiSC). Promoters that get methylated 
between serum ESC and EpiSC states are enriched in being bivalent. Lower panel shows collated CpG density and chromatin state of all 
promoters that become significantly differentially (hydroxy)methylated between at least one pairwise PSCs comparison. 

(C) ProfiLes of 5mC and 5hmC at Stelta and Tell in PSCs. Shaded area indicates promoter region. 

(D and E) Enrichment dynamics of 5mC (blue) and 5hmC (red) by (h)meDIP-seq across ground-state ESCs in 2i, metastable ESCs in serum, 
and primed EpiSCs at (D) distinct genomic landmarks or (E) gene network promoters. 
Error bars represent SD of target enrichment(s). CGI, biologically defined CpG island. 



"primed" genes continue to acquire 5hmC in EpiSCs and 
remain depleted of 5mC, which may reflect the require- 
ment to keep these genes competent for expression in 
EpiSCs. Finally, germline-speciflc genes become robustly 
methylated in EpiSCs, which correlates with a dramatic 
loss of promoter 5hmC and may "lock in" a silenced state 
(Figure 2E). These data indicate a dynamic interplay be- 
tween promoter 5mC and 5hmC that may contribute to 
the maintenance of each pluripotent state. 



Imprint Stability in PSCs 

We next asked how 2i conditions affect genomic imprints, 
which are generally resistant to hypomethylation in vivo 
and only become demethylated in PGCs. We found that 
enrichment of 5mC at imprints was highly comparable be- 
tween ESCs in 2i and ESCs in serum, suggesting that they 
are resistant to the global hypomethylation in ground-state 
cells (Figure 3A) (Leitch et al., 2013). Moreover, imprints 
were fully maintained in EpiSCs, with the slight increase 
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Figure 3. 5mC and 5hmC Enrichment at Genomic Imprints and Enhancers 

(A) Box plots show enrichment of 5mC and 5hmC over gametic-imprinted DMRs (n = 21) in PSCs. 

(B) Distribution of 5mC and 5hmC over the imprinted PeglO DMR (gray shading) in PSCs and Ell. 5 PGCs. 

(C) Bisulfite sequencing reveals loss of imprints at Peg3 and PeglO DMRs in independent earLy-passage XX ESC lines but imprint main- 
tenance in control EpiSCs and XY ESCs. 

(D) Clustering of PSCs according to 5mC and 5hmC enrichment at imprinted DMRs. EGCs and XX ESCs (ESC-F) cluster separately from other 
PSCs due to erasure of imprints. 

(E) Average enrichment of 5hmC at all gametic-imprinted DMRs in PSCs and Ell. 5 PGCs. 5hmC is enriched specifically in PGCs, where 
imprints are in the process of dynamic erasure. 

(F) The profiles of 5mC (upper) and 5hmC (Lower) over poised enhancers in EGCs and ESCs in 2i, male and female ESCs in serum, and EpiSCs 

(G) GluMS-qPCR of quantitative (%) 5mC and 5hmC levels at the Foxol, Gck, Klf4, and Nanog enhancers in ESCs (XY) cultured in serum or 2i. 
Error bars represent SD of biological duplicate experiments. **p < 0.01. 



in 5mC enrichment possibly a consequence of denser DNA 
methylation within each methylated allele in EpiSCs rather 
than an increase in allelic methylation. In contrast, we 
observed that 5mC was significantly depleted from many 
imprints in EGCs in 2i, which reflects their origin from 
PGCs (Figures 3 A and 3B). We also found imprints to be 
almost completely erased in our low-passage female ESCs 
(Figures 3A and 3B), which we confirmed in two further in- 



dependent XX ESC lines (Figure 3C). Thus, genomic im- 
prints are resistant to the globally hypomethylated state 
in 2i conditions, but like EGCs, XX ESCs are highly 
disposed to imprint erasure. Indeed, EGCs and XX ESCs 
cluster together according to 5mC enrichment of "im- 
prints," whereas EpiSCs and XY ESCs in 2i or serum cluster 
separately (Figure 3D). Because our XX ESCs exhibited no 
significant loss of promoter 5mC and only a modest 
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reduction in global 5mC, they apparently become hypo- 
methylated preferentially at genomic features that are nor- 
mally most resistant to DNA demethylation in vivo — IAP 
and genomic imprints — through an as yet undetermined 
mechanism. Notably imprints were highly depleted of 
5hmC in all cell types, which contrasts with PGCs, which 
exhibit highly significant 5hmC enrichment coincident 
with imprint erasure (Figure 3E) (Hackett et al., 2013). 
Interestingly, EGCs exhibited the greatest enrichment of 
5hmC at imprints despite reduced levels of 5mC substrate, 
which could be a consequence of their derivation from 
PGCs (Figure 3A). However, the low absolute levels of 
5hmC at genomic imprints in PSCs, in general (Figures 3 
and 2D), likely contribute to their relative resistance to 
demethylation in most contexts, including in the globally 
hypomethylated state of 2i conditions. 

Differential DNA Modifications at Enhancers 
We next tested the enrichment of 5mC and 5hmC at en- 
hancers in PSCs. Poised enhancers, denned by DNase 
I-hypersensitive sites (DNase HSs) and H3K4mel enrich- 
ment in ESCs, were depleted of 5mC but enriched for 
5hmC in ESCs in serum, as expected (Figure 3F) (Yu et al., 
2012). In contrast, active enhancers (additionally marked 
by H3K27ac) exhibited no enrichment for 5hmC and 
remained depleted of 5mC, suggesting that once activated, 
enhancers dispense with the 5hmC signature (Figure S3). A 
similar profile at enhancers was observed for XX ESCs, 
further underscoring the general similarity between XX 
and XY ESC methylomes on feeders (excepting imprints 
and IAP) (Figure 3F). Interestingly, in contrast to 5hmC 
enrichment in serum conditions, poised enhancers ex- 
hibited a significant depletion of 5hmC in both ESCs and 
EGCs in 2i. Moreover, whereas EpiSCs maintained the 
5hmC signature at this enhancer set, it was also associated 
with acquisition of 5mC (Figure 3F). These data suggest 
that cells in each pluripotent state may utilize or are pre- 
marked for the use of distinct enhancer sets, which may 
have roles in cellular competence and epigenetic priming 
(Song et al., 2013; Shen et al., 2013). 

To further investigate distal regulatory elements, we used 
glucosyltransferase methylation-sensitive qPCR (GluMS- 
qPCR) to examine poised enhancers linked with the Pro- 
B cell- and liver-specific genes Foxol and Gck, respectively, 
and active enhancers linked to the pluripotency genes Klf4 
and Nanog (Whyte et al., 2013). We found 5mC levels to be 
between 10% and 50% at all these enhancers in ESCs in 
serum, consistent with them being low-methylated regions 
(LMRs) that are associated with regulatory elements (Fig- 
ure 3G) (Stadler et al., 2011). At all enhancers examined, 
5mC was reduced 2- to 4-fold in 2i, which is proportion- 
ately in line with the global hypomethylation in ground- 
state cells. Moreover, consistent with hmeDIP-seq, the 



poised Foxol and Gck enhancers were highly enriched for 
5hmC (>5%) in ESCs in serum, whereas the active Klf4 
and Nanog enhancers were relatively depleted (<1%). How- 
ever, the enrichment of 5hmC at poised enhancers was 
dramatically lost in 2i (>10-fold), to a significantly greater 
extent than global 5hmC depletion (Figure 3G). Thus, 
5hmC is apparently a specific marker of poised lineage- 
specific enhancers in ESCs in serum, but not 2i. The enrich- 
ment of 5hmC may be important for resisting the high 
de novo methyltransferase activity in serum conditions, 
thereby antagonizing accumulation of repressive DNA 
methylation at enhancers that need to remain competent 
for downstream commitment. 

Molecular Events during Transition to Ground-State 
Pluripotency 

Because of the dramatic DNA hypomethylation in 2i condi- 
tions, we next investigated the kinetics of transition from 
serum to 2i. We initially shifted ESCs in serum to 2i condi- 
tions and tracked gene expression changes. As expected, 
naive pluripotency markers such as Prdml4, Esrrb, and 
Stella were upregulated upon transfer to 2i, whereas 
primed/early differentiation markers including FgfS, Lefty 1, 
and Bmp4 were strongly downregulated (Figure 4A) (Marks 
et al., 2012). We also observed that the de novo methyl- 
ation machinery including Dnmt3a, Dnm3b, and Dnmt3L 
was rapidly repressed in 2i, whereas maintenance genes 
such as Dnmtl and Uhrfl exhibited no change or were 
upregulated, respectively. Expression levels of Tetl were 
reduced by ~2-fold in 2i but remained expressed at very 
high overall levels, whereas Fet2 and Fet3 were not signifi- 
cantly affected, albeit the latter was expressed at very low 
absolute levels. Notably, changes in gene expression 
occurred rapidly upon transfer to 2i, with most changes 
being comparable to ESCs maintained in 2i (more than 
five passages) within 48 hr (Figure 4A). 

The rapid transition of a serum ESC gene expression 
profile to a profile comparable to 2i conditions prompted 
us to examine the kinetics of DNA methylation changes. 
We initially used bisulfite sequencing and found that 
LINE-1 and the promoters at Dazl, Rhox6/9, and ElfS 
showed some loss of DNA modification at 48 hr with signif- 
icant demethylation by 120 hr (Figure 4B). In particular, 
Rhox6/9 and Elf5, which are representative of relatively 
CpG-dense and CpG-poor promoters, respectively, were 
highly hypomethylated after 120 hr in 2i. In contrast, we 
observed no significant changes at IAP and the imprinted 
Kcnqlotl locus, consistent with our observation that these 
regions are resistant to demethylation in 2i. We next used 
GluMS-qPCR to track quantitative changes of 5mC and 
5hmC at a high temporal resolution at genes with CCGG 
motifs. We found that Dazl, Rhox6/9, and Aszl promoters 
underwent a slight reduction of 5mC within 24 hr with a 
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Figure 4. Molecular Dynamics during Transition to Ground-State Pluripotency 

(A) Gene expression Levels in ESCs in serum (FCS) and at 48 and 120 hr and more than five passages after switch to 2i conditions. Shown are 
absolute mRNA Levels as a proportion of Gapdh. 

(B) Bisulfite sequencing shows DNA methylation of ESCs in serum and at 48 and 120 hr after switch to 2i at gene promoters and LINE-1 
repeats. 

(C) Kinetics of quantitative (%) 5mC and 5hmC Levels upon transition of metastable ESCs in serum (FCS) to ground-state pluripotency (2i) 
assayed by GLuMS-qPCR. 

*p < 0.05 and **p < 0.01, paired Student's t test. Error bars are SD of independent duplicate experiments performed in technical triplicates. 
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Figure 5. PRDM14 Can Drive Partial DNA Demethylation 

(A) Induction of exogenous Prdml4 (+D0X) in ESCs in serum has significant effects of gene expression. Shown are expression Levels at 48 
and 120 hr relative to uninduced (-D0X) ESCs. 

(B) Western blot validates PRDM14 overexpression after DOX treatment and downregulation of DNMT3A and DNMT3B. 

(C) Quantitative 5mC and 5hmC Levels by GluMS-qPCR during induction of exogenous Prdml4 in ESCs in serum. 

(D) Bisulfite sequencing of Dazl and £7/5 from ESCs in serum that have uninduced Prdml4 (—DOX), induced Prdml4 for more than five 
passages (+D0X), or control uninduced in 2i conditions. 

*p < 0.05 and **p < 0.01, paired Student's t test. Error bars represent SD of biological duplicate experiments performed in technical 
triplicates. 



progressive loss thereafter. In parallel, the quantitative 
levels of 5hmC remained at least stable over the first 
72 hr, despite the reduction in its obligate 5mC substrate 
(Figure 4C). Indeed, we observed that 5hmC levels may 
actually increase modestly at these loci, with a peak be- 
tween 48 and 72 hr after shift to 2i conditions. Because 
there is no known 5hmC maintenance mechanism (Ficz 
et al., 2011) and ESCs transferred to 2i replicate normally, 
the observation that 5hmC is initially maintained, and 
possibly enriched, whereas 5mC is lost suggests that 
5hmC is being continually targeted to these loci at an 
increasing rate to account for the decreasing level of the 
underlying 5mC substrate. Indeed, plotting the ratio of 
5hmC:5mC at tested loci demonstrates the relative enrich- 
ment of 5hmC at these loci during transition to ground- 
state pluripotency (Figure S4). These data suggest that 
5hmC turnover and retargeting may play a role in medi- 
ating demethylation upon transfer to 2i. 

Prdml4 Contributes to DNA Demethylation in 
Ground-State ESCs 

To further investigate the mechanism that leads to global 
hypomethylation in 2i, we thus considered two possibil- 



ities, which may operate synergistically. First, that 5hmC 
is targeted to the genome at an increased rate during tran- 
sition to ground-state pluripotency, perhaps owing to 
changes in chromatin accessibility or enhanced TET 
activity, leading to faster turnover of 5mC and an equilib- 
rium shift toward hypomethylation. Second, the upregula- 
tion of Prdml4 observed in 2i (Figure 4 A) leads to a rapid 
repression of the de novo methyltransferase machinery, 
thereby promoting the passive loss of 5mC owing to ineffi- 
cient maintenance (Chen et al., 2003). Indeed, we and 
others have recently reported that Prdml4~ / ~ ESCs exhibit 
higher steady-state DNA methylation at tested loci than 
wild-type (Grabole et al., 2013; Yamaji et al., 2013). 

To examine whether Prdml4 directly mediates the dy- 
namic DNA demethylation effects of transfer to 2i, we 
generated ESCs carrying DOX-inducible expression of 
Prdml4. Induction of Prdml4 (+DOX) in serum resulted 
in downregulation of primed/early differentiation markers, 
likely due to repression of FGF signaling (Yamaji et al., 
2013; Grabole et al., 2013), and some upregulation of naive 
pluripotency markers by 48 hr (Figure 5A). Importantly, 
PRDM14 targets Dnmt3a, Dnmtib, and Dnmt3L all 
exhibited significant repression upon DOX treatment, 
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Figure 6. TET-Mediated 5hmC Contributes to DNA Demethylation in 2i 

(A) Gene expression dynamics of wild-type and Tetl/Tet2 knockout (DKO) ESCs upon transition to ground-state pLuripotency normalized to 
wild-type and DKO in serum. 

(B) Kinetics of 5mC and 5hmCin wild-type and DKO ESCs at 48 and 120 hrand more than five passages after shift from serum to 2i conditions. 
Shown are CpG-dense promoters (Dazl and Rhox6/9), CpG-poor promoters (Capn9 and Esrrb), introns {Cull), and repeats (LINE-1). 

(C) Bisulfite sequencing of Dazl, Elf5, and Rhox9 promoters and LINE-1 in wild-type and DKO ESCs after five passages in 2i conditions. 
*p < 0.05, Student's t test. Error bars are SD of independent duplicate experiments performed in technical triplicates. 



suggesting that the de novo methylation machinery is 
impaired upon Prdml4 expression, even in serum condi- 
tions (Figures 5A and 5B). Nonetheless, GluMS-qPCR 
revealed that whereas there was some DNA demethylation 
at Dazl, Aszl, and Rhox6/9 upon Prdml4 overexpression, 
the rate and extent were significantly lower than seen 
upon transfer to 2i conditions (compare Figure 5C with Fig- 
ure 4C); this despite ~ 20-fold Prdml4 upregulation 
compared with ~6-fold in 2i. This observation was 
confirmed by bisulfite sequencing (Figure S4). To test 
whether Prdml4 could promote demethylation over 
extended time periods, we induced Prdml4 for more than 
five passages. Under these conditions, we did observe 
significant demethylation relative to uninduced (— DOX) 
at tested loci; however, DNA methylation still remained 
at significantly higher levels than in ESCs shifted to 2i (Fig- 
ure 5D). This difference between Prdml4 overexpression 
and 2i may in part reflect the fact that Dnmt3a, Dnmt3b, 
and Dnmt3L are only partially repressed by Prdml4 over- 
expression alone but may also be indicative of other 



Prdml 4-independent DNA demethylation pathways being 
operative during transition to ground-state pluripotency. 
Because Prdml4~ / ~ ESCs retain DNA methylation in 2i 
(Leitch et al., 2013; Grabole et al, 2013; Yamaji et al., 
2013), but overexpression of PRDM14 does not fully 
recapitulate global hypomethylation, our results suggest 
that Prdml4 is necessary but not sufficient to induce a fully 
hypomethylated epigenome in ground-state pluripotent 
cells. 

TET-Mediated 5hmC Contributes to Global 
Hypomethylation 

We next turned to the role of 5hmC in promoting DNA 
demethylation by using Tetl and Tet2 compound double- 
knockout (DKO) ESCs (Dawlaty et al., 2013). Expression 
analysis suggested that upon shift to 2i, naive markers 
such as Prdml 4 and Esrrb showed no significant change 
in the rate or extent of upregulation in DKO ESCs relative 
to wild-type (Figure 6A). Moreover, Lefty 1 and Fgf5 were 
both rapidly repressed in control and mutant cells. 
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However, we did observe that methylation-dependent 
genes such as Dazl and Aszl were upregulated faster in con- 
trol cells as compared to DKO, implying an inhibition of 
DNA demethylation at these loci in DKO ESCs (Figure 6A). 

To further investigate this, we examined the DNA deme- 
thylation dynamics in DKO ESCs. In serum conditions, we 
observed no significant difference in DNA methylation 
levels between control and DKO ESCs by bisulfite 
sequencing. However, upon transfer to 2i, the rate of deme- 
thylation was impaired specifically in DKO cells (Figure S4). 
GluMS-qPCR analysis showed that CpG-dense promoters 
(Dazl and Rhox6/9), CpG-poor promoters (Esrrb), and exons 
(Cull) all exhibited a delay in the rate of 5mC depletion in 
DKO cells, whereas some loci (Capn9) were apparently un- 
affected (Figure 6B). As expected, 5hmC was undetectable 
in DKO ESCs at all loci, but a peak of quantitative 5hmC 
enrichment occurred in wild-type controls at 48 or 120 hr 
after shift to 2i. The relative delay of 5mC erasure in DKO 
cells was greater than the reported slight proliferation 
defect (Dawlaty et al., 2013), suggesting that these results 
do not only reflect slower DNA replication in DKO ESCs. 
Indeed, importantly, even after more than five passages 
in 2i conditions, DKO cells still exhibited significantly 
elevated levels of 5mC at multiple loci, suggesting that in 
the absence of TET1 and TET2, ESCs are intrinsically 
blocked from reaching complete global hypomethylation 
(Figure 6B). Bisulfite sequencing confirmed that Dazl, 
Elf 5, and Rhox6/9 promoters and LINE-1 repeats are all 
substantially more methylated in DKO cells relative to 
wild-type in 2i, albeit these loci have undergone significant 
demethylation relative to serum conditions (Figure 6C). 
Collectively, these results support an important role for 
TET-mediated 5hmC conversion in enhancing complete 
DNA demethylation in ground-state pluripotency. This 
likely occurs synergistically with repression of the de 
novo methyltransferases, in part by PRDM14, to promote 
a robust hypomethylated state. 

DISCUSSION 

Our findings document that discrete pluripotent states — 
naive ground state (2i), metastable (serum), and primed — 
are each associated with distinct 5mC and 5hmC profiles. 
Because these states are interchangeable, culture condi- 
tions thus dictate the epigenomic landscape of pluripotent 
cells. Notably, culture in 2i directs a globally hypomethy- 
lated state, with our genome-wide analysis indicating 
that most loci are susceptible to DNA demethylation, as 
recently suggested by Ficz et al. (2013) and Habibi et al. 
(2013). Exceptions are imprints and IAPs, which in accor- 
dance with their resistance to demethylation during preim- 
plantation reprogramming, remain methylated in 2i. We 



found that in general, 5mC is acquired during progression 
from ground-state pluripotency to primed pluripotency. 
However, whereas the levels of global 5mC dramatically in- 
crease between ESCs in 2i and those in serum, promoters 
are relatively resistant to acquisition of 5mC until progres- 
sion to the primed state, characterized by EpiSCs. In 
contrast, 5hmC is particularly enriched at promoters in 
ESCs cultured in serum, which may function to protect 
promoter regions from the high global 5mC and de novo 
methylase activity in this condition, thus maintaining a 
permissive transcriptional state. Consistently, this observa- 
tion is particularly prominent at pluripotency-associated 
genes. Moreover, poised, but not active, enhancers are 
highly enriched for 5hmC specifically in ESCs in serum, 
which may function to protect these regions from stable 
DNA methylation in the absence of functional activation. 
This could potentially maintain their competence to 
enable unbiased lineage commitment. Indeed, ESCs 
depleted of TET proteins exhibit skewed lineage differenti- 
ation (Koh et al., 2011; Dawlaty et al., 2013). 

The dual presence of high de novo methylase activity 
and high 5hmC activity in ESCs in serum may represent 
capture of a unique epigenetic state not present in vivo, 
whereby the methylome is reflective of a more advanced 
state of differentiation (hypermethylated) but retains 
competence for pluripotency through demethylation of 
key regulatory elements. Such targeted hypomethylation 
(at some enhancers and promoters) may explain the 
conundrum that ESCs in serum are associated with high 
levels of global DNA methylation normally linked with 
lineage restriction but still exhibit the pluripotent proper- 
ties of the hypomethylated preimplantation epiblast. In 
contrast, the striking global hypomethylation in 2i culture 
parallels the naive epiblast in vivo and supports the 
paradigm that 2i conditions more faithfully recapitulate 
the unrestricted state of the inner cell mass from which 
ESCs are derived. 

Notably, ESCs in serum can rapidly initiate DNA deme- 
thylation upon shift to 2i. Mechanistically, this is founded 
upon Prdml4 upregulation driving repression of Dnmt3 
genes and TET-mediated conversion to 5hmC. These 
processes are likely complementary and shift the equilib- 
rium between de novo methylation activity and 5hmC- 
mediated demethylation toward hypomethylation. Thus, 
upon transition to 2i, Dnmt3 enzymes are repressed, in 
part by Prdml4, which suppresses the potential for de 
novo DNA methylation activity at the same time as 
5hmC activity is at least maintained. This promotes greater 
demethylation both independent of 5hmC, through ineffi- 
cient maintenance methylation, and via 5hmC, owing to 
the absence of cyclical remethylation following 5hmC 
conversion. Indeed, it is possible that TET catalytic activity 
is actually enhanced in the 2i system either through altered 
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chromatin states, such as reduced H3K27me3/H3K9me3 
(Marks et al., 2012), or by a direct effect of culture condi- 
tions (Blaschke et al., 2013) because we observed a modest 
peak of 5hmC at 48-72 hr after transfer to 2i. In contrast, 
the reciprocal transition toward a primed state is associated 
with increased de novo activity through upregulation of 
Dnmt3 genes and parallel downregulation of Tetl (Fig- 
ure S4), which shifts the equilibrium toward stable 5mC 
maintenance. 

Knockout of Dnmta3 and Dnmt3b in ESCs leads to pro- 
gressive hypomethylation, supporting the notion that de 
novo activity is required to maintain 5mC in ESCs, possibly 
to offset persistent 5hmC conversion (Chen et al., 2003). 
However, importantly, Dnmt3a/b null ESCs only exhibited 
an ~10%-20% reduction in global 5mC after five passages 
compared with >2-fold in 2i (Jackson et al., 2004; Chen 
et al., 2003; Leitch et al., 2013), suggesting that ground- 
state conditions must activate additional DNA demethy- 
lation processes independent of Dnmt3a/b repression. 
Indeed, Tetl/Tet2 null ESCs used here exhibit a partial block 
to DNA demethylation in 2i conditions, supporting an 
important role for 5hmC conversion. Nonetheless, a signif- 
icant degree of hypomethylation was still observed in DKO 
ESCs, implying that DNA demethylation is not reliant on 
5hmC per se. Thus, a dual model may operate whereby 
repression of de novo methylation activity is sufficient to 
drive hypomethylation at slow kinetics, whereas the pres- 
ence of abundant 5hmC activity may enhance both the 
rate and the extent of demethylation, particularly at 
CpG-dense regions that are preferential TET-binding sites, 
such as the Dazl promoter (Williams et al., 2012). Notably, 
though, base-resolution data reveal that 5hmC is detect- 
able to some extent throughout the ESC genome (Yu 
et al., 2012), and thus, 5hmC may play at least a subordi- 
nate role in promoting hypomethylation genome wide. 
Indeed, we observed a partial suppression of demethyla- 
tion at introns (Cull) and CpG-poor promoters {ElfS, Essrb, 
and Capn9) in DKO ESCs in addition to CpG-dense regions 
{Dazl and Rhox9). This model of synergistic demethylation 
mechanisms is reminiscent of reprogramming phases in 
PGCs and the preimplantation embryo, where significant 
global 5hmC conversion has been reported, and that coin- 
cides with periods of suppressed maintenance and/or de 
novo methylation activity (Hackett et al., 2013; Kagiwada 
et al, 2013; Iqbal et al, 2011; Pastor et al., 2013). 

The functional consequences of global hypomethylation 
in 2i conditions are difficult to disentangle from the direct 
effects of ground-state pluripotency. However, one possible 
outcome is a relaxation on epigenetic barriers imposed by 
DNA methylation. For example, hypermethylation of the 
ElfS promoter has been reported to act as an epigenetic 
barrier to prevent ESCs in serum from entering into extra- 
embryonic lineages (Ng et al., 2008). Because we observed 



rapid and dramatic DNA demethylation and activation of 
ElfS expression after switching ESCs to 2i (Figure S4), this 
should enable these ESCs to contribute to extraembryonic 
tissues and thus be functionally "totipotent." Indeed, we 
found that ESCs carrying a reporter for constitutive- Venus 
expression and cultured in 2i robustly contributed to extra- 
embryonic tissues and that this property occurred soon 
after the switch in culture conditions (unpublished data), 
a finding supported by a recent study (Morgani et al., 
2013). Thus, global hypomethylation in 2i may contribute 
to generating a less-restricted in vitro state, which is closer 
to functional totipotency, possibly through demethylation 
of ElfS. Additionally, DNA hypomethylation per se appears 
to promote homogeneity among sister cells after ESC 
division and may therefore directly promote self-renewal 
(Jasnos et al., 2013). 

In summary, we have established the 5mC and 5hmC 
profiles of distinct PSCs and revealed that DNA demethy- 
lation during transition to ground-state pluripotency 
is directed by synergistic TET-mediated 5hmC and 
PRDM14-directed repression of Dnmt3a/Dnmt3b/Dnmt3L. 
The globally hypomethylated state of cells in 2i may be 
an in vitro correlate to the preimplantation epiblast or 
migrating PGCs, and thus, transfer to 2i conditions may 
provide a tractable system for mechanistic studies into 
developmental epigenetic transitions. 



EXPERIMENTAL PROCEDURES 

Cell Lines and Culture Conditions 

Animal studies were authorized by a UK Home Office Project 
License and carried out in a Home Office-designated facility. Cell 
lines used for (h)meDIP-seq profiling were derived from embryos 
generated from crossing mixed background Ocf4APE-GFP trans- 
genic males with 129SvEv strain females and cultured as follows. 
For 2i culture, cells were maintained on laminin (10 ng/ml) in 
N2B27 basal media supplemented with 1 uM PD0325901 (MEK in- 
hibitor), 3 nM CHIR99021 (GSK3 inhibitor), and LIF at 10 ug/ml 
(prepared in house). For serum conditions, cells were maintained 
on a MEF feeder layer in DMEM-F12 (GIBCO) supplemented 
with 15% FCS, 0.1 mM MEM nonessential amino acids, 2 mM 
L-glutamine, 1 mM sodium pyruvate, 0.1 mM |3-mercaptoethanol, 
and LIF at 10 ng/ml. Cells were passaged by dissociation with 
trypsin and replating every 2-3 days. Prior to (h)meDIP-seq 
analysis, cells were purified from feeders by sequential preplating 
for 2 x 45 min. For primed conditions, EpiSCs were maintained 
on fibronectin in DMEM-F12 with N2 and B27 supplement 
(Invitrogen) media with 2 mM L-glutamine, 0.1 mM MEM nones- 
sential amino acids, 0.1 mM (3-mercaptoethanol, 50 ng/ml BSA, 
supplemented with recombinant human activin A (20 ng/ml; 
PeproFech), and bFGF (12 ng/ml; Invitrogen). EpiSCs were 
passaged and replated using Accutase (PAA) every 1-2 days. For 
serum to 2i transfer experiments, feeder-free ESCs were derived 
on gelatin in GMEM supplemented with 15% FCS, 0.1 mM MEM 
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nonessential amino acids, 2 mM L-glutamine, 1 mM sodium pyru- 
vate, 0,1 mM p-mercaptoethanol, 100 U/ml penicillin/100 ng/ml 
streptomycin, and LIF at 10 ng/ml and shifted to N2B27 supple- 
mented with 1 uM PD0325901, 3 uM CHIR99021, 1% knockout 
serum replacement (KSR), and 10 ng/ml LIF on fibronectin. Upon 
switching to 2i, cells were passaged every 2-3 days by dissociation 
with TrypLE with no obvious effects on cell viability or the ex- 
pected proliferation rate. 

(h)meDIP-Seq, Bisulfite Sequencing, and GluMS-qPCR 

Genome-wide (h)meDIP-seq, bisulfite sequencing, and GluMS- 
qPCR were performed as described (Hackett et al., 2013). Full 
experimental procedures are detailed in the Supplemental Experi- 
mental Procedures. 
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